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ABSTRACT 



O I The radial component of the peculiar velocities of galaxies cause 

Q I displacements in their positions in redshift space. We study the effect of the 

^ I peculiar velocities on the linear redshift space two point correlation function. 

I Our analysis takes into account the radial nature of the redshift space distortions 



and it highlights the limitations of the plane parallel approximation. We 



I consider the problem of determining the value of /? and the real space two 

^ I point correlation function from the linear redshift space two point correlation 

^ I function. The inversion method proposed here takes into account the radial 

nature of the redshift space distortions and can be applied to magnitude limited 
redshift surveys that have only partial sky coverage. 



Q ■ Subject headings: Galaxies: Clustering - Large Scale Structure of the Universe 

methods: analytical 



^ ' 1. Introduction. 



The redshift of a galaxy has information about both its position and the radial 
component of its peculiar velocity, and as a consequence galaxy-galaxy correlations in 
redshift space differ from the correlations amongst their real positions. Using linear theory 
and the plane-parallel approximation. Kaiser (1987) showed that the redshift space power 
spectrum is the real space power spectrum amplified by the factor (1 + fig^yU^)^ where /i is 
the cosine of the angle between the line of sight and the wave vector. He also pointed out 
that the anisotropy of the redshift space power spectrum can be used to measure the value 
of Qq. 



^Present address: Department of Physics, IIT Kharagpur, 721 302, India, e-mail: 
somnath@phy. iitkgp . ernet . in 
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The plane-parallel approximation assumes that the pairs of galaxies between which the 
correlation is being measured are sufficiently far away so that their separation subtends 
a very small angle at the observer and the displacements in redshift space caused by 
their peculiar motions may be treated as parallel. Most of the subsequent work (Lilje & 
Efstathiou 1989; McGill 1990; Loveday et al. 1992; Hamilton 1992; Gramman, Cen & 
Bahcall 1993; Bromley 1994; Fry & Gaztaiiaga 1994; Fisher et al 1994; Fisher 1995; Cole, 
Fisher & Weinberg 1994, 1995; Matsubara & Suto 1996;) are based on the plane-parallel 
approximation. 

A proper analysis however requires that the radial nature of the displacements in 
redshift space is taken into account. This is required if pairs of galaxies which subtends 
a large angle are to be also included in the analysis. This was first taken into account in 
the analysis of Fisher, Scharf & Lahav (1994) who decomposed the angular behaviour of 
the density field (in redshift space) into spherical harmonics and integrated out the radial 
dependence after multiplying the density field with a Gaussian radial window function. The 
power spectrum of the coefficients of this expansion was then used to determine the value 
of j3 = ^Q^/fe, and they obtained the value (3 = 1.0 ± 0.3 on applying this technique on 
the 1.2-Jy /i?^ 5 redshift survey. This method was improved by Heavens & Taylor (1995) 
who expanded the radial dependence of the density field into spherical Bessel functions, 
but both these analysis have the drawback that they require a prior knowledge of the linear 
power spectrum P{k). These methods were further refined by Ballinger, Heavens & Taylor 
(1995) who do not fix the shape of P{k) but allow it to vary in six bins in k space. All these 
methods have the limitation that they require the galaxy survey to have full sky coverage 
and they cannot be applied to two dimensional redshift surveys. 

Hamilton & Culhane (1996), and Zaroubi & Hoffman (1996) have calculated the 
linear two point correlation in redshift space taking into account the radial nature of the 
distortions. Zaroubi Sz Hoffman (1996) have also investigated the mode-mode couphng that 
arises when the analysis is done in Fourier space, but they have not addressed the problem 
of determining the value of j3 and the real space correlation in detail. This problem has been 
discussed in detail by Hamilton & Culhane (1996) who propose a method for determining f3 
in a manner which does not require any prior assumption about the real space correlation 
function. They have studied the eigenfunctions of the 'spherical distortion operator' which 
relates the real space correlation to its redshift space counterpart. The observed redshift 
space correlations is expanded in terms of these eigenfunctions and the ratio of the values 
of these coefficients can be used to obtain /3. The problem of dealing with the effect of the 
radial nature of the redshift space distortions on the two point correlation function has also 
been considered in a recent paper by Szalay, Matsubara & Landy (1998). 



-3- 



In this paper we have investigated the effect of the radial nature of the redshift space 
distortions on the linear two point correlation function . The redshift space two point 
correlation is a function of the triangle formed by the observer and the pair of galaxies for 
which the correlation is being measured. In order to get a better understanding of the effects 
of the redshift space distortions we study in some detail how the redshift space correlation 
function changes with the shape of the triangle, the value of /3 and the slope of the real space 
correlation. We also address the question as to when the radial nature of the distortions is 
important and when they may be ignored and the plane parallel approximation be used 
instead. In section 3. we address the problem of determining the value of (5 and the real 
space correlation function from the observed redshift space correlation function taking into 
account the radial nature of the distortions. The analysis of Hamilton & Culhane (1996) is 
based on the assumption that the selection function has a power law spatial dependence. 
We have investigated whether this assumption is justified for a magnitude limited sample. 
In this paper we discuss the inverse problem for two different situations (1.) assuming that 
the selection function is a power law (2.) for a more realistic form of the selection function 
which can be used in a magnitude limited survey. Both the inversion methods proposed 
here can be applied to redshift surveys that have partial sky coverage. 



2. The linecir redshift space correlation. 

A large amount of our information about the spatial distribution of galaxies is inferred 
from redshift surveys which provide angular positions and redshifts of a large number of 
galaxies. The distance to galaxies is very hard to measure and the analysis of redshift 
surveys has to rely on the redshift as an indicator of the distance to the galaxies. This has a 
drawback because in addition to the Hubble expansion the redshift has contributions from 
the radial component of the peculiar velocity of our Galaxy and the galaxy being observed. 
The peculiar velocity of our Galaxy has been determined from the dipole anisotropy 
observed in the CMBR (Kogut et al. 1993) and this can be corrected for in all observations 
made from our Galaxy but the contribution from the radial component of the peculiar 
velocity of the galaxy being observed remains in its redshift. 

Using the vector z to denote the position of a galaxy in the 3-dimensional redshift 
space formed by the angular positions and redshifts, the relation between z and the actual 
position of the galaxy x can be written as 

z = x + x(v-x). (1) 
Here the hat denotes a unit vector (x = x/x) and v • x is the radial component of the 
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peculiar velocity of the galaxy, The units have been chosen so that units where the speed of 
light c, and the present value of the scale factor a and its time derivative a are all unity, 

The problem is how to use quantities measured from the distribution of galaxies in 
redshift space to draw inferences about the actual distribution of the galaxies. In addressing 
this problem it is also necessary to take into account the fact that usually the galaxies in a 
redshift survey are not selected uniformly from the region of space surveyed. For magnitude 
limited surveys the selection criteria is a function of the actual distance from the observer 
and it is represented by the selection function 0(x) which gives the fraction of galaxies 
selected in the survey as a function of the distance from the observer. 

Taking these effects into account, the observed number density of galaxies in 
redshift space n^{z) can be related to the number density of galaxies in real space 
n(x.) =< n > (1 + 5(x)) and the peculiar velocity field v(x) as (Kaiser 1987) 



n^(z) = $(^) <n> {1 + 5{z) 



z \ d \og{z) J 



z . v(z) . (2) 



where 9^ = z ■ Vz is used to denote the derivative in the radial direction. 

This equation is valid at linear order in v, and it has the assumption that v << z and 
is a slowly varying function. It is convenient to define a function 

"«^(-^) 

where a{z) = 2 for a volume limited sample for which the selection function is a constant. 

In the linear regime, in the presence of only the growing mode of perturbations 
(Peebles 1980), it is possible express the perturbation and the peculiar velocity in terms of 
a potential as 

5(x) = VV(x) and v(x) = -/3V^/'(x) . (4) 

Here (3 = Q^'^/b where Qq is the density parameter and b is the bias parameter which takes 
into account the fact that the galaxies may be a biased tracer of the underlying matter 
density which determines the peculiar velocities. 

We use these relations in equation (^ to write the number density of galaxies at the 
point z in redshift space as 



$(2) < n > < 1 + 



a[z\ 



V'(z) 



(5) 
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and we use this to calculate the linear redshift space two point correlation function 
^^{7,1,7.2) defined as 

where the angular brackets < ... > denote ensemble average. In evaluating this we encounter 
the ensemble average < ■?/'( 2:1 )■?/'( 2:2) >. As the universe is statistically homogeneous and 
isotropic we can define 0(-22i) =< i'{'^i)i'{'^2) > which is a function of only the magnitude 
of the vector Z21 = Z2 — zi. Using this we obtain 



(^2i). (7) 



for the linear redshift space two point correlation function. 

This should be compared with the real space two point correlation function 

e(zi, Z2) =< 5(zi)5(z2) >= VV(^2i) . (8) 

which depends only on ^21- 

Equations (|^) and (j^) together are equivalent to the expression for the linear redshift 
space correlation function derived by Hamilton and Culhane (1996). Equation (H) can be 
inverted to relate various derivatives of the potential which appear in equation (|^) to 
integrals of ^(x) and this is described in Appendix A. 

The expression for the redshift space two point correlation function presented here 
is valid in the regime where ^(zi,Z2) <^ 1. In addition there are the restrictions that the 
redshifts Zi and Z2 are in a range where they are much larger than -\/< f ^ > (the r.m.s. 
peculiar velocity) and where the selection function does not vary too rapidly. 

Unlike the real space two point correlation ^(| zi — Z2 |) which depends on just on 
the distance between the points zi and Z2, the redshift space counterpart depends on the 
triangle formed by the observer O and the points zi and Z2, and we next investigate this 
behaviour in some detail. 

The behaviour of ^^(zi,Z2) is relatively simple in the situation where the two edges 
of the triangle zi and Z2 are made very large keeping Z21 fixed . In this limit z'l and z'2 
are nearly parallel i.e. lim^j^oo z"i = z"2 = n and the peculiar velocities of the galaxies at 
Zi and Z2 can be treated as being parallel. In addition, if the selection function is such 
that \miz^oo{.c>i{z) / z) = 0, then the terms involving a(z) can be dropped and equation (0) 
becomes 

e^(zi,Z2)= [V^^,+/3(n- Vz^J^l' 0(^21) (9) 
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and we have the hnear redshift space two point correlation in the plane-parallel 
approximation (PPA). In this limit the redshift space two point correlation depends on the 
length of just one edge of the triangle Z21 = Z2 — zi, and it depends on on the angle between 
Z21 and the line of sight n. This angular dependence introduces anisotropy in the redshift 
space correlation and this is well understood in PPA (Hamilton 1992). Here we investigate 
the behaviour of the redshift space correlation function in a more general situation where 
the plane-parallel approximation cannot be applied and the radial nature of the distortions 
has to be taken into account. In the rest of the discussion in this section we use the value 
a = 2 which corresponds to a volume limited sample where the selection function is a 
constant. 

In order to separately study the dependence of ^^(zi, Z2) on the shape and the size of 
the triangle formed by the observer O and the points zi, Z2 we consider a situation where 
the real space correlation function has a power law behaviour ^{x) oc x~'^ . In this case the 
effect of changing the size of the triangle is very simple C,^{yzi,yz2) = y^''C,^{zi,Z2), and 
the ratio 

W(Z1,Z2)=^^(Z1,Z2)/^(Z21) (10) 

depends only on the shape of the triangle. We have used the function w(zi,Z2) to study 
how ^^(zi, Z2) varies with the shape of the triangle. 

We parameterize triangles of all possible shapes by first carrying out the following 
operations which leave ty(zi,Z2) unchanged: (1) Label the larger of the two sides that 
originate from O as zi (2) Rotate the triangles around O so that they all lie in the x-y 
plane with zi along the x axis (3) Scale the triangles so that zi = 1. At the end of these 
operations, for all the triangles zi corresponds to a unit vector in the x direction while Z2 
lies in the xy plane and it is restricted to be inside a circle of unit radius centered around O. 
Triangles which lie in the lower half plane can be related to triangles with the same shape 
in the upper half plane by reflecting Z2 on the x axis. It is thus possible to parameterize 
triangles of all possible shapes by using the vector Z2 which is restricted to lie inside the 
upper half of a circle of unit radius centered around O. Figure |I]. shows the observer O at 
the point (0,0), the point zi at (1,0) and the semi-circle shows the region inside which the 
point Z2 must lie. Every point in the semi-circle corresponds to a triangle with a different 
shape, and one possible triangles is shown in figure 0. We have used this parameterization 
to study how w{z2) varies with the shape of the triangle and the results of this study are 
presented in the form of contour plots which show contours of equal w plotted at equal 
intervals of w for triangles of all possible shapes. 



We have studied the behaviour of the function w{z2) for three different cases 7 = 4, 3.5 
and 2.5 with (3 = 1 and the corresponding contour plots are shown in figures |2^, ^ and 
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respectively. The value of w{z2) along the 45° cuts shown in figures 2a, 3a and Ja are 



plotted in figures |2E, and HB, and the purpose of these graphs is to show the value of w 



corresponding to the different contours. Figures ^ and ^ also show w{z2) for other 
values of (3 for which we have not shown contour plots. 

Going back to the contour plots, we expect PPA to be valid in the region near the 
lower right hand corner of the figures where Z2 is nearly equal to Zi and Z21 ^ Zi. For 
comparison we have calculated w{z2) using PPA for the case with 7 = 4 and /3 = 1 for 
triangles of all shapes (i.e. also beyond the region where we would expect PPA to be valid) 
and this is shown in figure 2.c. A comparison of figures 2. a and 2.c shows that as expected 
the figures match in the region around Zi, but we start seeing considerable differences from 
the predictions of PPA as either (1). the angle between zi and Z2 is increased or (2). the 
sizes of zi and Z2 start to differ significantly. We find that PPA correctly describes the effect 
of the redshift space distortions on the linear two point correlation function provided the 
two vectors zi and Z2 do not differ by more than 30 % i.e. 2:21 < .3 Zi. The function w{z2) 
behaves quite differently from the predictions of PPA once the two redshift space vectors zi 
and Z2 differ by more than 30 % and the radial nature of the distortions become important 
in this situation. 

A comparison of figures 2. a, 3. a. and 4. a shows how the effect of the redshift space 
distortion changes with the slope of the real space correlation 7. For 7 = 4 and 3.5 the 
redshift space correlation is of the same sign as the real space correlation in the forward 
direction (when zi and Z2 are in the same direction) and the redshift space correlation 
function changes sign as Z2 is turned away from zi, whereas for 7 = 2.5 the behaviour is 
just the opposite. We also find that the shapes of the contour lines changes significantly as 
the value of 7 is changed. For all the cases the effect of redshift space distortions becomes 
stronger as Z2 approaches the observer and w{z2) diverges in the limit Z2 ~^ 0. This 
behaviour is due to the factor a{z2)/z2 which appears in equation (|^, and the cause of the 
divergence can be related to the fact that under the map from real space to redshift space a 
non-zero volume element in real space collapsed to a point at 2; = 0. 

In figure 2.d we have shown w{z2) for the case where 7 = 4 and (3 = 0.2. Comparing 
this with figure 2. a we see that a change in the value of (3 can make a qualitative difference 
in the behaviour of the redshift space two point correlation function. Although the shape 
of the contour lines is similar in both the figures, the region where w{z2) is negative differs 
significantly in the two plots. 
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3. Inverse Problem. 

In this section we address the problem of using the linear redshift space correlation 
function measured from redshift surveys to determine the value of j3 and the real space 
correlation function. 

The inverse problem is very easily solved if ,^^(zi,Z2) is restricted to a region where 
PPA is valid. The function (^^(zi, Z2) can then be written as a function of Z21 =| Z2 — zi | 
and yU which is the cosine of the angle (shown in figure |1| ) between Z21 and fi - the line of 
sight to the pair of galaxies. In PPA the angular dependence of ^^{z2i,fi) is very simple 
and it can be expressed in terms of Legendre polynomials as 

00 

e^(z21,/i) = E^K-")e/^(^2l) (11) 
1=0 

where the ^P'{z2i)s are different angular moments of ^^(^21, A*)- As shown by Hamilton 
(1992), only the first three even moments have nonzero values and the value of (3 can be 
obtained from the ratios of these angular moments determined from redshift surveys. 

The situation is changed if the radial nature of the redshift space distortion is taken 
into account as the behaviour of ^^(zi, Z2) now depends on the shape of the triangle formed 
by zi,Z2 and O. This has been studied by Hamilton & Culhane (1996) who have expressed 
the shape dependence of ^^(zi, Z2) in terms of five shape functions and they have proposed 
a method for measuring the value of (3 based on this. In their work Hamilton and Culhane 
(1996) have made the simplifying assumption that the selection function can be described 
by a power law ^{z) oc which implies that a{z) = a is a constant. 

Here we propose a different method for analyzing the linear redshift space two point 
correlation function and determining the value of (3. In the first part of the analysis 
presented below we assume that a{z) = a is a constant and later on in a separate subsection 
we treat the inverse problem for a more realistic selection function. 

In this analysis zi and Z2, the two sides of the triangle which originate from O, are 
not treated on a equal footing. Writing one of the sides (say zi) as zi = zih, all possible 
triangles can be parameterized using the lengths of two of the sides - zi and Z21, and /i 
the cosine of the angle (shown in figure |l|) between Z21 and fi. This way of parameterizing 
is very similar to that used in PPA, except that in PPA fi is the common direction along 
which both Zi and Z2 lie whereas now fi refers to the direction of one of the sides Zi. Also, 
in PPA two parameters Z21 and /i suffice to describe the behaviour of whereas we 
require an additional third parameter zi if the radial nature of the distortions are taken 
into account. For a fixed value of 2:21, PPA corresponds to the limit Zi —* 00, and hence the 
value of Zi does not appear in in PPA. 



- 9 - 



It is convenient to express zi in terms of the dimensionless ratio 

s - — 



(12) 



and we use (s,2;2i,/i) to parameterize For a fixed value of the hmit s — > 
corresponds to PPA and the effects of the radial nature of the redshift space distortions 
become more important as s increases. 

The set of parameters (5,^21, /u) has the feature that any triangle can be described by 
two sets of values of s and /i. For example for a triangle formed by a pair of galaxies (which 
we call A and B) and O, we get different values of s and yU depending on whether we label A 
as Zl or we label B as Zi. Here we take the point of view that both the labellings should be 
used and as a consequence each pair of galaxies contributes to ^'^s, Z2i,ft) for two different 
values of s and /i. 

We first consider the angular dependence of ^^(5,221,/^), and following the analysis 
used in PPA (equation |lT]) we decompose in terms of Legendre polynomials. We find 
that unlike in PPA, now all the moments have non-zero values and this does not provide a 
convenient way of analyzing The analysis becomes considerably simpler if we use the 
redshift weighted correlation function C,^ defined as 



P(S, Z21, /i) = ^ C"(S, Z21, fx) . 



(13) 



where {z2/zi)'^ can be written as (1 + + 2s/i). The angular dependence of ^^{s, Z21, fx) is 
much simpler, and expanding it in terms of Legendre polynomials 



P(s,Z21,/^) = Ee^(^'-22l)Pn(/i) 



(14) 



n=0 



we find that only the first five moments have non-zero values which are shown below 

2 ^ 1 



^o{s,Z2i] 



ii{s,Z2i] 



ms,Z2i) 



l + -f3+-[3')az2i) + s 



(1 + ^/?+ 1/5^)^(^21) 



-a'(3%{z2i) + -f3{-6 + a-2(3-a(3 + a'(3)Uz2i) 

D 9 



16 



6 



28 



2 + + -13') e(^2i) - [-f^ + -f3') 6(^21 

+ ■^^"/^ -^a/5^i(^2i) - ^6(^21) - 7^/5(3 - a)^4{z2i] 
lb 6 15 



4 4 



(e(^2l)-6(^2l))+S^ 



-m + p)az2i) 



(15) 



(16) 
(17) 
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+ ^f3 (-6 -ia + 2p + ai3- a^p) ^2{z2i) - (4 + 
iiis,Z2i) = s^p[{l + P)^{z2i)-{l-PM^2i)-2PU^2i)_ 



(18) 



(19) 



We have used the relations given in the Appendix A to express the derivatives of the 
potential ^(-221) in terms of volume integrals of the real space two point correlation function 
which are defined as 



The first point to note is that for a fixed value of 2:21, in the limit s ^ we recover the 
results calculated by Hamilton (1992) in the plane parallel approximation. As expected, 
the odd moments all vanish in this limit. 

The effect of including the radial nature of the distortions manifests itself through the 
s dependent terms. The different angular moments all have a very simple dependence on 
s and this involves at most a cubic polynomial in s. The angular moment has no s 
dependence, and the expression calculated for this moment using PPA remains unchanged 
if the radial nature of the distortions is taken into account. 

The procedure for determining the value of f3 from redshift surveys is now quite 
straightforward in principle. The first step is to estimate the redshift weighted correlation 
function ^^{s, Z21, ij) using all the pairs of galaxies in the survey. For a fixed value of 2:21, 
for a survey which extends from a redshift Za to a redshift z^, the variable s will lie in the 
range Z2i/za >s> Z2i/zb. 

The second step is to decompose the angular dependence of the observationally 
determined ^^(s, ^2i,A*) in terms of Legendre polynomial. This is possible only for those 
values of (s, Z21) where there are observations of ^^{s, Z21, /j) for both positive and negative 
II, and in general the range of (s, Z21) where this is possible will depend on the geometry of 
the redshift survey. The angular decomposition can be used to obtain estimates for the first 
five angular moments ^^(s, 2:21) to ^4^(5, 2:21) ■ We expect all the higher moments to be zero 
in the linear regime. 

For a fixed value of 2^21, the s dependence of the angular moments ^^(s, 2^21) is very 




(20) 



and for n > 1 




(21) 
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simple in the linear regime. The third step is to do a least squares fit for the s dependence 
of the Cni^y -2^2i)s using polynomials in s of the form predicted by equations (|T5D to (pJ|). For 
example, the monopole determined from redshift surveys can be fitted using a quadratic 
function of the form ^1^(5,2:21) = ^(^(^21) + ^02(^21) s^, where C,qq{z2i) and are the 
unknown quantities which have to be determined from the fit. 

The coefficient of the terms s° in the fits for the even moments gives ^(f (0, Z21), ■^2^(0, -221) 
and ^^{s, Z21) which correspond to the values of these moments in PPA. Effectively, this 
allows us to use pairs of galaxies for which the radial nature of the redshift space distortion 
is important (large values of s) to predict the behaviour of the angular moments at s = 
which corresponds to PPA. Once these are known we can use the method proposed by 
Hamilton (1992) for determining /? and ^(-221) using results calculated in the plane-parallel 
approximation. 

There are various possible ways of doing this and one way is to use the relations 



:i + lp+\p')k{z2i)-u^2i) 

3 5 L 

(^/?+^/?2)[e(^2i) -6(^21) 



eo>,^2i)-4/ ^oiO,rydr (22) 
Z21 Jo 



e?(o,^2i). 



(23) 



The right hand side of equations ( |22D and (^) involve quantities which can be determined 
from the redshift surveys by using the procedure discussed above. The ratio of these two 
equations can then be used to calculate j3, which can in turn be used to determine ^{x). 



3.1. A Realistic Selection Function. 

Until now we have discussed the inverse problem for a situation where the selection 
function has a specific form ^{z) oc which implies that a{z) = a is a constant. In this 
subsection we consider a realistic situation where we have a magnitude limited survey. We 
first investigate the validity of the assumption that a{z) is a constant. 

For a magnitude limited sample with a lower apparent magnitude limit rrimin and an 
upper magnitude limit rrimax, the selection function is related to N{M) the differential 
galaxy luminosity function by 

/•Mmaxiz) 

<^{z) = Ci N{M)dM (24) 



where the limits in the integral can be related to the limits in the apparent magnitude using 
the relation between absolute and apparent magnitudes, and the galaxy luminosity function 
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N{M) is usually modeled using the Schechter function (Schechter 1976). The normalization 
constant Ci in equation (|2^ ) is of no interest in this discussion as it does not affect a{z), 
and it is only the shape of the selection function which is of interest. 

As an example we have used the apparent magnitude limits and the luminosity 
function given for the N112 subsample of the Las Campanas Redshift Survey - LCRS 
(Lin et. al. 1996) to calculate ^{z) and a{z) (shown in figures ^ and P) for go = -5 and 
-f^o = lOOkm/s/Mpc . We have set the normalization Ci to an arbitrary number chosen for 
the convenience of plotting ^{z). 

We find that the assumption that a{z) is a constant does not correctly describe 
the behaviour of the calculated values of a{z) shown in figure @. The behaviour of 
a{z) in the redshift range 0.05 < z < 0.19 can be fitted using a function of the form 
a{z) = a — z^ /q^ with a = 1.7 and q = 0.072. This corresponds to a selection function of 
the form ^{z) = Cz"~^e~^^/^'^^ (shown in s figure P) and we find that this fits the calculated 
selection function to better than 2% in the range 0.05 < z < 0.19. Figure (|^) also shows 
n{< z) - the fraction of the galaxies we expect to find at redshifts less than z, and we see 
that about 90% of the galaxies are expected to lie in the redshift range over which this fit 
is valid. 

In the rest of this paper we assume that the selection function is of the form 

$(z) = C2"-2e-"'/2g2 (25) 

with 

a{z) =a- z^/q^ (26) 

where a and q are constants which have to be determined for the particular galaxy survey 
being considered. We expect the selection function to behave as predicted by equation ( p5|) 
for in any magnitude limited survey where the luminosity function is of the Schechter form 
which is a product of a power law and an exponential. It may also be noted that a fit 
of the form a{z) = —z^/q^ which corresponds to a selection function ^{z) = Ce~^^^'^'^^ is 
reasonably accurate (~ 10%) and this may also be used instead of equation (|25|). 

The form of the selection function in equation (123) introduces a new length-scale q and 



^'^(zi,Z2) now behaves quite differently from the case where a{z) is a constant. The most 
important difference is that now the limit limz^oo / z diverges and hence ^'^(zi,Z2) as 
given by equation (|^) diverges in the plane parallel approximation which corresponds to the 
limit zi — > oo and Z2 — > oo with Z21 held fixed. Before proceeding further it is necessary to 
clarify the fact that this divergence in the behaviour of ^^(zi, Z2) is not a physical effect and 
we do not expect to find a divergent behaviour in the ^^(zi,Z2) determined from redshift 
surveys when Zi or Z2 are made very large. The divergent behaviour predicted by equation 
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(0) arises because the derivation of this equation involves Taylor expanding $(z + t>p) {vp 
is the radial component of the peculiar velocity) in powers of vp. For the selection function 
given in equation this involves expanding the function e'^^+^^^^Z^"?^ = Q-{z'^+'2zvp+vj,)/2q^ 



in powers of vp. Such an expansion is valid only if z vp -C q"^ and vp <^ q, and the expansion 
is invalid if we take the limit z —>■ oo. As a consequence equation (0) also becomes invalid 
for very large values of zi or Z2, and this gives rise to a divergent behaviour in ^'^(zi,Z2). 
Equation is valid provided y/ < v"^ >Z\lq^ <^ 1, \/ < v"^ >Z2/q^ <^ 1 and ^(2:21) is in 
the linear regime, and it does not correctly describe the behaviour of if any of these 
conditions are not satisfied. 

We next check over what redshift range equation (|^) is valid for LCRS. Using the 
values V< f ^ > ~ 0.004 and q = 0.07, we find that equation (j^) can be used at linear 
scales over the entire redshift range for which the fit given by equation (^) is valid. Having 
ascertained the fact that equations and (p5|) are both valid over the redshift range in 
which most of the galaxies lie, we proceed to investigate the inverse problem for a situation 
where the linear is described by these equations. 

As a consequence of the divergent behaviour discussed above, it is not possible to 
extrapolate the determined from pairs of galaxies for which equations (|^) and are 
valid to obtain in the plain parallel approximation. Hence, it is not possible to obtain 
the value of P using the method discussed earlier for the case where a{z) is a constant, and 
a different inversion method has to be used for the selection function given by equation 



We find that in this case it is convenient to use (zi, -221,/^) to parameterize the redshift 
weighted two point correlation function 



^2 

^""{z^, Z21, fi) = ^^""{z,, Z21, fi) (27) 
^1 



where Z2 can be written as Z2 = zf + Z21 + 2ziZ2i^- 

Decomposing C,^{zi, Z21, fi) in terms of Legendre Polynomials (equation |14D we again 



find that only the first five moments are non-zero. In this case it is more convenient to 
work directly in terms of the potential 4>{z2i) instead of expressing the derivatives of 0(2:21) 
in terms of volume averages of ■^(^21)- The expressions for the angular moments and 
are given below 



e?(^l,^2l) 





(28) 
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35 



dz- 



2 Z21 



21 



dzo 



■(P{z2l) 



(29) 



21 



The expressions for the other non-zero ^^s are very similar to equation ( p8| ) and these are 
presented in Appendix B. The equation for ^^{zi,Z2i) is much simpler than the equations 
for the other angular moments and it does not involve Zi. This equation can be integrated 
once to obtain d^(j){z2i)/dz2i in terms of the ^^(^21) determined from the redshift survey, 
and this gives us the relation 

' ^ ^^e^-Z-'V" rfye-^V^^ef(l/) + e(^--^i)/^^^ (30) 

) Jzr. 



dzii 



^21 



where we only use the values of ^^^(2:21) at separations Z21 > zi where we expect the 
correlations to be in the linear regime. The constant of integration Ai corresponds to 
d'^(p{z2i)/dz2i at the point 2:21 = zl. Equation (|30D can be integrated once more to obtain 

(iV. . 35 



dz^ 



21 



/?2 



Jzl J Z]^ Jzl 



(31) 



where A2 is a constant of integration which corresponds to the value of (I){z2i) / dz2i at 
Z21 = Zl- Equation ( pO| ) can also be used to obtain the relation 



''ef(.2O + ^^0(^2O. 



(32) 



21 



These equations (^), (0) and (|3^) allow us to determine the value of the derivatives 
of the potential (f){z2i) in terms of the observed ^f^(-22i) and three unknown parameters (3, 
Ai and A2. Using these relations in equation (^Sf) we obtain the expression for ^|^(-2i, ^21) 
presented below 



(^1, Z21] 



+ 



+ 



+ 



+ 



7^21 

2zi 

7/1 

^1 \(3 

7Z2l ( 

2g2 1^ 

5 



1 



ef(^2i) 



+ 



^21 



4g2 13 q 



a 

Zl 



Zl_ 



Hdre''/"' r dy 

Jzl Jzl 



21 

5zi 



a 

X 



Zl 



Ai 



Z21 



4+-(4 + /?)U2e(^--^')/'?' 



(33) 
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This equation allows us to use the ^^^(-221) determined from observations to predict the 
values of ^^(21,^21)- The quantities which involve the observed values of ^^{zi,Z2i) have 
been enclosed in square brackets in the above equation. The relation between the observed 

and the predicted values for also involves three parameters (3, Ai and A2 which 
are the only unknown quantities in equation (^). By comparing the values of C,3{zi, Z21) 
predicted by equation (|33|) with the observed values of C,^{zi, Z21) it is possible to determine 
the value of these parameters /3, Ai and A2 for which the predicted best fits the observed 

A similar procedure can also be carried out using the angular moments ^q, and 
and the expressions for these quantities in terms of the observed C,^ and the parameters (3, 
Ai and A2 are presented in appendix B. The best fitting values of the parameters f3, Ai 
and A2 determined using the four different angular moments should be consistent and this 
provides a way of checking the validity of the method proposed here. 

One the values of the parameters (3, Ai and A2 are known it is quite straight forward 
to determine the real space correlation C(^2i) = ^^0(^21)- 

4. Summary and Discussion. 

We have studied how the peculiar velocities distort the linear two point correlation 
function in redshift space. Our analysis takes into account the radial nature of the 
distortion. We have compared this with the linear two point correlation calculated in 
the plane parallel approximation and we find that there are significant differences in the 
behaviour of ^^(zi,Z2) if the two redshift space vectors zi and Z2 differ by more than 
30 %. The effect of the radial nature of the redshift space distortions become important 
when either the angle between zi and Z2 becomes large, or the lengths of zi and Z2 differ 
significantly and the plane parallel approximation does not correctly describe the effect of 
redshift space distortions under these curcumstances. 

We have also addressed the problem of extracting the value of /? and the real space 
correlation function from the redshift space correlation function taking into account the 
radial nature of the distortions. This problem has been studied earlier by Hamilton & 
Culhane (1996) who have assumed that a{z) is a constant. We have tested the validity of 
this assumption for a magnitude limited survey and we find that such an assumption is not 
justified. The inversion scheme proposed by Hamilton & Culhane (1996) can be applied 
only to volume limited samples where a{z) =2. 

In the first part of our analysis we have followed Hamilton & Culhane (1996) in 
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assuming that a{z) is a constant. The inversion procedure proposed here is quite different 
from that proposed by Hamilton & Culhane (1996). The main difference is that we 
propose the use of a redshift weighted two point correlation function instead of the 
correlation function The function has the advantage that if we decompose its 
angular dependence in terms of Legendre polynomials only the first five angular moments 
are non-zero, and all the higher angular moments are zero. This is not the case if we use 
instead. The expressions we obtain for the angular moments of are very closely related 
to the results one gets in the plane parallel approximation and this makes the proposed 
inversion procedure very simple. 

For a magnitude limited sample the selection function can be very well approximated 
by ^{z) — C^r^-^e"^^/^*^ and we have analyzed the inverse problem for such a situation. 
We find that in this situation also the redshift weighted correlation function has only five 
non-zero angular moments, and the procedure proposed for determining /3 in section 3.1 is 
based on this. This procedure can be applied to magnitude limited samples. It also has the 
feature that it uses the values of the redshift space correlation only at separations larger 
than some separation Zl which can be chosen so that all scales larger than z^ are in the 
linear regime. This has the advantage that the inversion procedure uses only the values of 
the redshift space correlation function from scales which are definitely in the linear regime 
and the inversion procedure is not affected in anyway by the non-linear scales. 

Both the inversion schemes discussed in this paper can be applied to redshift surveys 
with partial sky coverage - for example they can be applied to the Las Campanas Redshift 
Survey where the observations are restricted to six thin conical slices. 



A. Appendix A. 

The equation 

= V'(t>{z) (Al) 

can be inverted to obtain 

- fay)y'dy. (A2) 

D Z Jo 

where Zi refers to the components of z, di refers to the components of Vz and C is a 
constant whose value is fixed by the boundary conditions. It is most natural to choose the 
condition 

lim didj(t){z) = (A3) 
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which gives us 

C = - i{y)ydy. (A4) 

^0 



We also have 



1 /'^ 1 1 /"-^ 

didjdk(p{z) = --didjdkiz) i{y)y'^dy - -^didjdk{-) i{y)y*dy (A5) 



and 

1. ... ,1 



didjdkdi(t){z) = ^^^^^^{z) - -didjdkdi{z) f ^{y)y^dy 

I Jo 

- Ud^dkdii-) r ay)y'dy (A6) 

6 z Jo 



B. Appendix B. 



Here we present the expressions for the angular moments C,o{zi,Z2i), C,i{zi,Z2i) and 
^2{zi,Z2i) of the redshift weighted linear two point correlation function ^^{zi, Z21, fi) 
calculated for the selection function described by equation (pS]). 
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Equations (pOl 
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and (|32D can be written as 
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(B4) 
(B5) 
(B6) 



where the functions f2{z2i) and /3(2:2i) are defined as 

/3(2;2i) 

/2(^2l) 



2i 



and ^2(^21) and ^3(^21) are integrals of the observed (-221) 



^^2(2:21) 





rz2i 
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(B8) 
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Using equations (p^ ), (p5|) and ( pBD for the derivatives of 0(-22i) in equations (pT| ), 
( P2|) and ( |B3| ) , we obtain the following expressions for ^q, C,i and in terms of the 
observed and three parameters P, Ai and ^2. 
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Fig. 1. — This figure shows the parameter space corresponding to triangles of all possible 
shapes. The observer O at the origin and the point Zi at (1,0) form a fixed edge of the 
triangles. The third point Z2 can lie anywhere in the semi-circle, and different values of Z2 
corresponds to triangles of different shapes. One possible triangle is shown in the figure. 
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Fig. 2a. — This shows contours of equal w plotted in the parameter space corresponding to 
triangles of all possible shapes. The solid lines show positive values of w and the dashed 
lines show negative values of w. 




Fig. 2b. — This shows w along a cross-section at 45° through figure ^ plotted as a function 
of r -the distance of the point Z2 from the observer. Positive values of r correspond to 
the situation where Z2 and zi are in the same direction, and r is negative when they lie in 
opposite directions. 
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Fig. 2c. — This shows contours of equal w calculated using the plane parallel approximation 
plotted in the parameter space corresponding to triangles of all possible shapes. The solid 
lines show positive values of w and the dashed lines show negative values of w. 
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Fig. 2d. — This shows contours of equal w plotted in the parameter space corresponding to 
triangles of all possible shapes. The sohd lines show positive values of w and the dashed 
lines show negative values of w. 
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Fig. 3a. — This shows contours of equal w plotted in the parameter space corresponding to 
triangles of all possible shapes. The solid lines show positive values of w and the dashed 
lines show negative values of w. 
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Fig. 3b. — This shows w along a cross-section at 45° through figure 3a 
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Fig. 4a. — This shows contours of equal w plotted in the parameter space corresponding to 
triangles of all possible shapes. The solid lines show positive values of w and the dashed 
lines show negative values of w. 




Fig. 4b. — This shows w along a cross-section at 45° through figure 4a 
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Fig. 5. — This solid curve shows the selection function $(z) for the N112 sample of LCRS. 
The dashed-dotted curve shows the fit using equation (|25|). The dashed line shows n{< z) - 
the fraction of galaxies that are expected to lie below a redshift z. 




